E(5) and X(5) shape phase transitions within a Skyrme Hartree-Fock + BCS approach 
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Self-consistent Skyrme Hartree-Fock plus BCS calculations are performed to generate potential 
energy curves (PEC) in various chains of Pd, Xe, Ba, Nd, Sm, Gd, and Dy isotopes. The evolution 
of shapes with the number of nucleons is studied in a search for signatures of E(5) and X(5) critical 
point symmetries. It is shown that the energy barriers in the PECs are determined to a large extent 
by the treatment of the pairing correlations. 
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I. INTRODUCTION 

The ground states of atomic nuclei are characterized 
by different equilibrium configurations which correspond 
to different geometrical shapes. The study of these equi- 
librium shapes, as well as the transition regions between 
them, has been the subject of a large number of theo- 
retical and experimental studies (for a review, see, for 
example, Ref. [l| and references therein). 

Within the framework of algebraic models, the differ- 
ent nuclear phase shapes are put in correspondence to dy- 
namic symmetries of some algebraic structure that links 
a specific mathematical symmetry with a specific nuclear 
shape. Dynamic symmetries provide a useful tool to de- 
scribe properties of different physical systems since they 
lead to exactly solvable problems and produce all results 
for observables in explicit analytic form. 

In nuclear physics, the algebraic structure of relevance, 
according to the Interacting Boson Model (fBM) 0], is 
given by U(6). There are three dynamical symmetries 
characterized by U(5), associated to spherical symme- 
try, SU(3), associated with axially deformed symmetry, 
and S0(6), describing 7-unstable shapes. Experimental 
examples of all three types of symmetries have been rec- 
ognized in many nuclei. 

The phase shape transitions correspond to the break- 
ing of these dynamic symmetries and they occur as the 
number of nucleons change in the nucleus. Understand- 
ing the behavior of systems undergoing a phase transi- 
tion is of especial relevance, since a complicated interplay 
of competing degrees of freedom occurs at the critical 
points. 

lachello [1, Q introduced the E(5) and X(5) critical 
point symmetries within the framework of the collective 
Bohr Hamiltonian |5| under some simplifying approxi- 
mations. Critical point symmetries provide parameter 
free (up to scale factors) predictions of excitation spec- 
tra and electric quadrupole B(E2) strengths for nuclei at 
the critical point of a phase/shape transition. The geo- 
metrical shape of the ground state can be described [1] 



'Electronic address: ' sarriguren@iem. of mac . csic . es | 



by three Euler angles and by the quadrupole deforma- 
tion parameters (3 and 7. At the critical point of the 
phase transition the potential in the [3 degree of freedom 
can be approximated by a simple square well potential, 
which is decoupled from the potential in the 7-variable. 
In the case of the E(5) critical point symmetry, which 
corresponds to the transition from spherical vibrational 
U(5) to deformed 7-unstable 0(6), the potential is flat 
in the 7-direction. In the case of the symmetry X(5), 
related to the transition from U(5) to axially symmetric 
prolate SU(3), a harmonic oscillator potential is used in 
the 7-direction. 

Empirical evidence of these transitional symmetries at 
the critical points were soon found in ^'^^Ba [6^ for E(5) 
and in ^^^Sm J] for X(5). Other nuclei have been also 
identified as good candidates for those symmetries. Such 
is the case of ^^"^Vd 0, and i^s-isoxe [9], which pro- 
vide examples of E(5) symmetry. Other N=90 isotones 
like i5°Nd [13, ^^^Gd [lH, and i^Soy [12] also provide 
examples of X(5) symmetry. 

Algebraic models are very suitable for systematic stud- 
ies because they provide powerful predictions with a very 
small number of parameters, but in order to deepen into 
the details, one has to perform microscopic investigations 
of shape transitions and critical points which are, to a 
large extent, still missing. In particular, it is interest- 
ing to examine whether the assumptions of relatively flat 
potentials in E(5) and X(5) are justified in different mi- 
croscopic models and the self-consistent mean-field ap- 
proximation, based on parametrizations widely used all 
over the nuclear chart, appears as a very attractive ini- 
tial tool for a link between algebraic models and micro- 
scopic theories. In this context, the relativistic mean-field 
framework has been employed, in calculations of PECs 
as functions of the quadrupole deformation [l^, [ij, [3 ■ 
These studies have been performed for isotopic chains 
in which the occurrence of critical point symmetries has 
been predicted. Since flat PECs are one of the expected 
characteristics of critical point symmetries, constrained 
calculations in those isotopic chains should result in rela- 
tively flat PECs for nuclei with the critical symmetry. It 
has been shown that particular isotopes exhibit relatively 
flat PECs over an extended range of the deformation pa- 
rameter. Nevertheless, the behavior of the PECs, and 
particularly the potential barriers, are quite sensitive to 
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the relativistic interaction used [13| , as well as to the pair- 
ing treatment. Therefore, the question arises whether a 
similar situation occurs when using nonrelativistic effec- 
tive interactions to study candidates for critical point 
symmetries. Such systematic nonrelativistic studies are 
still missing and it is very interesting to compare the con- 
clusions extracted from them with those obtained using 
both algebraic models and relativistic mean-field approx- 
imations. 

In this work, we extend the calculations mentioned 
above to the case of nonrelativistic self-consistent Skyrme 
Hartree-Fock -I- BCS mean-field calculations. The shape 
phase transitions corresponding to E(5) and X(5) sym- 
metries are investigated systematically in various isotopic 
chains containing some of the suggested critical nuclei. 
In particular we study Pd, Xe, and Ba isotopic chains as 
candidates to E(5) symmetry, and Nd, Sm, Gd, and Dy 
isotopic chains as examples of X(5). 

The paper is organized as follows. In Sec. lUwe present 
a brief description of the theoretical formalism (Hartree- 
Fock -I- BCS) used to obtain the main ingredient of the 
present study, i.e., the PECs for the considered isotopic 
chains. For a more detailed account the reader is referred 
to the corresponding literature. Sec. IIIII contains our 
results with a discussion on the sensitivity of the PECs to 
the effective nucleon-nucleon force and to the treatment 
of the pairing correlations. Sec. IIVI is devoted to the 
concluding remarks. 



II. THEORETICAL FRAMEWORK 

The microscopic approach used in this work (i.e., 
HF+BCS) is based on a self-consistent formalism built on 
a deformed Hartree-Fock (HE) mean-field, using Skyrme 
type energy density functionals. Pairing correlations be- 
tween like nucleons are included in the BCS approxi- 
mation. It is well known that the density-dependent 
HF+BCS approximation provides a very good descrip- 
tion of ground-state properties for both spherical and 
deformed nuclei [TBj and it is at present one of the state- 
of-the-art mean- field descriptions (l7j . 

There are two leading choices of methods to solve the 
deformed HF-I-BCS equations. One option is the use 
of a coordinate space mesh. In this case one solves 
the HF+BCS equations for Skyrme type functionals via 
discretization of the individual wave functions on a 3- 
dimensional cartesian mesh [l^ . This corresponds to an 
expansion on a specific basis of Lagrange polynomials 
associated with the selected mesh. The other common 
choice is to expand the single particle wave functions into 
an appropriate orthogonal basis (commonly the eigen- 
functions of an axially symmetric harmonic oscillator po- 
tential). In the present study we perform calculations 
with both methods. We use the code ev8 [3] in the first 
case and follow t he p rocedure based on the formalism 
developed in Ref. in the second case. 

In this work we consider the parametrization SLy4 (20| 
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FIG. 1: Potential energy curves obtained with SLy4 and fixed 
pairing gaps, using different methods for solving the HF+BCS 
equations. 
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FIG. 2; Potential energy curves obtained with three different 
Skyrme forces. For details see the main text. 



of the Skyrme force in the particle-hole channel, although 
we also show results in some instances for the forces Sk3 
[2l| and SG2 They are examples of global effec- 

tive interactions of Skyrme-type that have been designed 
to fit ground state properties of spherical nuclei and nu- 
clear matter properties. While SkS is the most simple 
one, involving in particular a linear dependence on the 
density, SLy4 is one of the most recent parametrizations 
of Skyrme forces. 

As we move away from closed shells, pairing correla- 
tions play an important role ^23.] and should be taken into 
account. If one were dealing with a fundamental many- 
body Hamiltonian, one would proceed to apply Hartree- 
Fock-Bogoliubov formalism to it. However, dealing with 
Skyrme forces that have been simplified with the aim of 
reproducing average or bulk properties of the nucleus, 
one would have to include additional parameters in or- 
der to guarantee that sensible pairing matrix elements 
are obtained [13, [Hi- 
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FIG. 3: Potential energy curves obtained with different treat- 
ments of pairing. For details see the main text. 



In this study pairing correlations are taken into ac- 
count in the BCS approximation. Several options have 
been investigated. Our main option will be the use of a 
zero-range density-dependent pairing force [23 |. 



y(ri,r2) = -g(l-P'^) 



l-^)^(r,-r.), (1) 



where P'^ is the spin exchange operator, p{r) is the nu- 
clear density, and pc — 0.16 fm^'^. The strength g of 
the pairing force [Eq.(IT])] is taken g = 1000 MeV-fm^ 
for both neutrons and protons and a smooth cut-off of 5 
MeV around the Fermi level has also been used. Let us 
mention that very recently the parametrization SLy4 has 
been successfully applied in combination with the pairing 
interaction [Eq.dl])] (with g — 1000 MeV-fm^) in system- 
atic studies of correlation energies from ^^O to the super- 
heavies [25,] and in global studies of spectroscopic prop- 
erties of the first 2+ states in even-even nuclei [26| . This 
is the main reason for selecting the combination SLy4 in 
the particle-hole channel and the interaction [Eq. ([1])] 
(with g = 1000 MeV-fm^) in the pairing channel as the 
leading choice for the present study. Additionally, results 
with the strength g — 1250 MeV-fm'^ will also be shown 
in some cases. 

Another common practice to include pairing correla- 
tions is to introduce a schematic seniority pairing force 
with a constant pairing strength G [1^ , we call this 
treatment constant-force approach, 



m,m'>0 



t t 



(2) 



The strength of the pairing force for protons and neu- 
trons Gp^n is chosen in such a way that the experimental 
pairing gaps extracted from binding energies in neigh- 
boring nuclei are reproduced. One can further simplify 
the pairing treatment by parametrizing the pairing gaps 



Ap.„ directly from experiment, we call this treatment 
constant-gap approach. The pairing strength and pair- 
ing gap are related through the gap equation [23j 



G'^^UuV^ 

y>0 



(3) 



where Vi, are the occupation amplitudes. 

The PECs shown in this study are computed micro- 
scopically by constrained HF+BCS calculations [l^, [H, 
[27| . These PECs are obtained by minimizing the corre- 
sponding energy functional under a quadratic constraint 
that holds the nuclear quadrupole moment fixed to a 
given value. In this work we show the energies as a func- 
tion of the (axially symmetric) mass quadrupole moment 
Q, which is related with the quadrupole deformation pa- 
rameter by the expression 



Q 



5 ^ < r2 > 



(4) 



in terms of the mean square radius of the mass distri- 
bution < >. Tables I and II contain the quadrupole 
parameters (32 corresponding to the equilibrium configu- 
rations for the isotopes considered in this work. 



TABLE I: Quadrupole deformation parameters (32 for the 
ground states in Pd, Xe, and Ba isotopes. SLy4 force and 
two values of the pairing strength g [MeV fm''] are used. 



Pd 


96 98 100 102 104 106 108 110 112 114 


g=1000 
g=1250 


0.10 0.12 0.15 0.17 0.18 0.17 0.17 0.18 0.19 -0.19 
0.00 0.10 0.14 0.15 0.16 0.16 0.15 0.13 0.13 0.11 


Xe 


118 120 122 124 126 128 130 132 134 136 


g=1000 
g=1250 


0.15 0.29 0.29 0.23 0.20 0.17 0.14 0.11 0.08 0.00 
0.09 0.27 0.27 0.20 0.18 0.15 0.13 0.09 -0.01 0.00 


Ba 


120 122 124 126 128 130 132 134 136 138 


g=1000 
g=1250 


0.34 0.28 0.27 0.25 0.23 0.20 0.17 0.13 0.11 0.00 
0.31 0.29 0.27 0.24 0.22 0.19 0.15 0.12 0.08 0.00 



III. DISCUSSION OF RESULTS 

It is known that PECs are sensitive to the effective 
nuclear force in both relativistic [31 and nonrelativis- 
tic I28L I29I approaches, as well as to pairing correla- 



tions 2^2^. Thus, it is worth starting our discussion 
on PECs by studying this sensitivity in our case. We 
perform this study on the example of ^^*Xe and ^^^Sm, 
which are proposed candidates for E(5) and X(5) sym- 
metries, respectively. 

In Fig. [l]we study the effect of using different methods 
of solving the HF+BCS equations, coordinate lattice [l^ 
and deformed harmonic oscillator basis [l^ . In these cal- 
culations we use the same SLy4 force in the particle-hole 
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TABLE II: Quadrupole deformation parameters /32 for the 
ground states in Nd, Sm, Gd, and Dy isotopes. SLy4 force 
and pairing strength g=1000 [MeV fm"^] are used. 



Nd 


142 144 146 148 150 152 154 156 




0.00 0.10 0.14 0.20 0.28 0.29 0.29 0.30 


Sm 


144 146 148 150 152 154 156 158 




0.00 0.07 0.16 0.20 0.25 0.29 0.30 0.31 


Gd 


146 148 150 152 154 156 




0.00 -0.05 0.17 0.22 0.25 0.28 


Dy 


148 150 152 154 156 158 




0.00 -0.08 -0.14 0.21 0.26 0.28 



channel and the same constant-gap treatment of pairing. 
The results indicate that the PECs obtained are practi- 
cally the same in these two cases and in general in all the 
isotopes we have considered in this work. 

In Fig. [2] we consider the PECs obtained from three 
different Skyrme forces, namely SLy4, SG2, and Sk3, us- 
ing the same pairing treatment. We can see that the 
location of the oblate and prolate minima appear at the 
same deformation no matter what the force is. How- 
ever, the relative energies of the minima and the energy 
barriers between the minima can change by a few MeV, 
depending on the force. 

In Fig. [3] we show the effect of pairing using SLy4 force 
in all cases. We consider the constant-gap approach and 
the constant-force approach. Again, we observe that the 
minima appear at about the same deformations but the 
energy barriers can change considerably depending on 
the approach used. 

From the analysis of Figs. [U [2] and [3] we conclude that, 
at least in the mass region studied, the PECs are not 
sensitive to the method employed to solve the HF+BCS 
equations (3D cartesian lattice or deformed harmonic os- 
cillator). We also conclude that the qualitative behavior 
of the energy profiles remains unchanged against changes 
in the Skyrme and pairing interactions in the sense that 
the deformations at which the minima occur are rather 
stable. Nevertheless, the relative energies of these min- 
ima, and particularly the energy barriers between them 
are very sensitive to the details of the calculation, espe- 
cially to pairing. This is also the case in relativistic mean- 
field calculations. The sensitivity of the PECs to the 
force can be seen in Ref. [13], where several forces were 
compared (NLl, NL3, NLSH, and TMl). The sensitivity 
to the pairing in relativistic calculations is also apparent 
if one compares the energy barriers in Sm isotopes ob- 
tained with the parametrization NL3 in Ref. [13], where 
a constant-gap approach was used, with those in Ref. 
p^ . where a pairing based on the Brink-Boeker part of 
the Gogny force was used. One should remark that, at 
variance with Ref. [l^ja self-consistent treatment of the 
pairing is used in Ref. [l5| . The main effect is a reduction 




Q[b] 

FIG. 4: Potential energy curves in Pd isotopes obtained from 
constrained HF-I-BCS calculations with the force SLy4 and a 
zero-range pairing force with two different strengths. 



of the energy barriers in the self-consistent treatment. In 
the case of nonrelativistic calculations with the Gogny in- 
teraction, since the same force is used in both mean-field 
and pairing channels, a fully self-consistent treatment is 
possible within a Hartree-Fock-Bogoliubov formalism. It 
will be very interesting to explore whether this reduction 
of the barriers is also present. Work along this line is in 
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FIG. 5: Same as in Fig. 4 but for Xe isotopes. 



progress. 

To appreciate the effect on the energy barriers in more 
detail, we show in Figs. [HISl and [6] the results for PECs 

96-ii4pj^ ii8-i36xe, and i^o-issgg^ isotopic chains, 
respectively. The quadrupole deformations correspond- 
ing to the ground states of these isotopes can be found in 
Table I. These chains contain E(5) candidates found from 
systematic studies on available data on energy levels, E2, 
El, and Ml strengths 0, i, i, [M, [Mi . 
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FIG. 6; Same as in Fig. 4 but for Ba isotopes. 



The PECs shown in Figs. 21 O and [5] have been com- 
puted within the mean-field scheme based on the Paris- 
Brussels code ev8, using the parametrization SLy4 and 
two different choices of the strength g of the zero-range 
pairing force [Eq.©], g = 1000 MeV-fm^ and g = 1250 
MeV-fm^ . From these figures one can see once more that 
oblate and prolate minima do not change significantly 
with the strength of pairing in each isotope. The only 
exceptions worth to mention are the cases ^^Pd, ^^"^Pd, 
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FIG. 7: Energy barriers in Xe isotopes corresponding to the 
cases shown in Fig. 5. 



and ^^''Xe, where we obtain the ground state at differ- 
ent deformations depending on the pairing strength g. 
This happens because of the flatness of the corresponding 
PECs, but in any case the energy difference between the 
minima are always very smafl. On the other hand, the po- 
tential barriers are clearly lower when pairing is stronger. 
We can also observe the transition from the spherical iso- 
topes 96pd (N=50), i36xe (N=82), and "Sfla (N=82) to 
the 7-unstable isotopes ^^^Pd, ^^^Xe, and ^^°Ba. We 
identify the isotopes los.nopd, i28,i30xe, and ^^^'^^"^Ba 
as transitional nuclei with rather flat PECs (especially 
with g=1250). These results confirm the assumed square 
well potential in the f} degree of freedom that leads to 
the critical point symmetry E(5). They are also in good 
agreement with the results obtained in Ref. ^] using 
the parametrization NL3 of the relativistic mean-field La- 
grangian and a pairing force based on the Brink-Boeker 
part of the Gogny interaction. 

Fig. [7] shows the evolution of the energy barriers as a 
function of the number of neutrons for Xe isotopes for the 
two values of the pairing strength g = 1000 MeV-fm'^ and 
g = 1250 MeV-fm^ in [Eq.((l|)]. As it can be observed, 
the barriers obtained with the strength g = 1250 MeV- 
fm^ are systematically lower than the ones obtained with 
g = 1000 MeV-fm3. 

Now we turn the discussion to several rare-earth iso- 
topic chains, Nd, Sm, Gd, and Dy, where some of the 
nuclei (N=90 isotones) have been identified as exhibiting 
X(5) behavior, a transition between a spherical shape 
and a well deformed prolate shape. For example, the nu- 
cleus ^^^Sm [?! was the first identified as exhibiting X(5) 
behavior and this is also the case for ^^'^Nd TO!]. Other 
candidates have also been sugg ested with further work 
(see, for example, Refs. [H, HI, H^l ) . 

Fig. [Hlshows the results for i42-i56]n^(J isotopes (Z=60). 
The corresponding ground state deformation parameters 
for these nuclei are given in Table [TTl We can observe 
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FIG. 8: Constrained HF-^BCS calculations in "2-156^^^^ 
topes with SLy4 and a zero-range pairing force with g = 1000 
MeV-fm"'. For a better comparison, the energies are shifted 
by 1 MeV for each isotope added, starting from A=142. Thick 
lines correspond to the isotopes that have been suggested to 
show a critical point symmetry. 
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FIG. 9; Same as in Fig. [8l but for Sm isotopes. 



a clear shape transition from spherical ^'^^Nd (N=82) to 
clearly prolate i^^-ise^d. The isotopes "^Nd (N=88) 
and ^^°Nd (N=90) show a transitional behavior with a 
rather flat minimum on the prolate side and additional 
minima on the oblate sector at 3 MeV and 4 MeV exci- 
tation energy, respectively. The energy barrier is about 
5 MeV in ^^SNd and 8 MeV in ^^ONd. 

Fig. [9] contains the PECs for ^'''^^i^^Sm isotopes 
(Z=62). The corresponding ground state deformation 
parameters are given in Table [TTl Similar arguments to 
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oblate and prolate minima in the case of Sm isotopes. 
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FIG. 11: Same as in Fig. |8] but for Gd isotopes. 



those in Nd are valid. The transitional behavior between 
spherical, ^^*Sm, and well prolate deformed, ^^'^^^^^Sm, 
appears again for the N=88 and N=90 isotopes, ^^°Sm 
and ^^^Sm, respectively. The energy barriers and energy 
differences between oblate and prolate minima in Sm iso- 
topes are shown in Fig. 1101 They are qualitatively similar 
to the barriers obtained in the other rare-earth isotopic 
chains. Our PECs for ^^''"^^^Sm agree qualitatively with 
those of Refs. and They also agree well with the 
results obtained in Ref. [35] with the Nilsson-Strutinsky 
+ BCS calculations. 

The same is true for Gd and Dy isotopes, whose PECs 
are shown in Figs. [11] and [T^ respectively. The ground 
state deformation parameters for the considered Gd and 
Dy isotopes are given in Table [TTl In this case, we find 
the transitional behavior for A=152,154 in Gd isotopes 
and A=154,156 for Dy isotopes. 



FIG. 12: Same as in Fig. |8l but for Dy isotopes. 

From the results described above we conclude that our 
calculations do not predict flat PECs for the N=90 iso- 
topes, in agreement with previous studies (see, for exam- 
ple, Ref. [il). 



IV. CONCLUSIONS 

In this paper we have investigated the potential energy 
curves in a series of isotopic chains containing transitional 
nuclei that exhibit a critical point symmetry behavior. A 
microscopic approach based on deformed HF-I-BCS cal- 
culations with Skyrme forces has been used. 

We have studied the sensitivity of our results to the 
effective nucleon-nucleon interaction and to the pairing 
force. From this analysis we conclude that the energy 
barriers in the PECs depend strongly on the details of 
calculations, especially on the pairing force. This is in- 
deed relevant for this work because our purpose is to 
study up to what extent the PECs exhibit a flat behav- 
ior at the critical point symmetries. 

We have found that the assumptions of flat potentials 
in the E(5) critical point symmetry is supported by the 
present microscopic calculations in los.iiop^j^ I28,i30^g^ 
and ^^^'^•^^Ba, that have been suggested as examples of 
E(5). In the case of X(5) we find that the rare-earth 
isotopes with N=88,90 show a transitional behavior that 
could be interpreted in terms of X(5) symmetry. How- 
ever, we do not find a flat behavior, in agreement with 
previous calculations [H, [3, . 

There remains a long list of tasks to be undertaken in 
the near future but the present study could be considered 
as a first step into a much more systematic exploration of 
the relation between algebraic models and (nonrelativis- 
tic) microscopic models. In particular, it is very impor- 
tant to understand how well the predictions of effective 
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interactions with predictive power all over the nuclear 
chart compare with those of the already mentioned alge- 
braic models. 
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